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Abstract 

We use the scale of Besov spaces B" T (0), a > 0, 1/t = a/d + 1/p, p fixed, 
to study the spatial regularity of the solutions of linear parabolic stochastic partial 
differential equations on bounded Lipschitz domains O C R d . The Besov smoothness 
determines the order of convergence that can be achieved by nonlinear approximation 
schemes. The proofs are based on a combination of weighted Sobolev estimates and 
characterizations of Besov spaces by wavelet expansions. 
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1 Introduction 

In this paper, the spatial Besov regularity of the solutions of linear stochastic evolution 
equations on bounded Lipschitz domains is studied. We combine regularity results by Kim 
[30] on stochastic partial differential equations (SPDEs, for short) on Lipschitz domains in 
terms of weighted Sobolev spaces with methods used in Dahlke, DeVore [13], where the 
Besov regularity of (deterministic) elliptic equations on Lipschitz domains is investigated. 
Our considerations are motivated by the question whether adaptive and other nonlinear 
approximation methods for the solutions of SPDEs on Lipschitz domains pay off in the 
sense that they yield better convergence rates than uniform methods. Thus referring to a 
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numerical theme and combining concepts and methods from different areas and scientific 
communities, the article is addressed to readers of both worlds: stochastic analysis and 
numerical analysis. Therefore, we give a rather detailed account in the first part of the 
paper, emphasizing conceptual and notational clarity. 

Our setting is as follows. On a finite interval [0, T] C [0, oo) let (w^) te r Q)T -i, k6M = 
{1,2,...}, be independent, one-dimensional standard Brownian motions with respect to a 
filtration (Ftjte&T] of cr-algebras on a complete probability space (fl, J 7 , P). Throughout 
the paper we assume that (FtjtefaT] is normal, i.e. the filtration satisfies the usual hypothe- 
ses, see, e.g. [18j Section 3.3.]. Let O C R d be a bounded Lipschitz domain. We consider 
the model equation 

d oo 

du= J2 aflUu ^ dt + J2d K dw?, u(0, • ) = M , (1) 

fl,l/ = l K=l 

for t G [0, T] and x G O. Here du is Ito's stochastic differential with respect to t, 
{a' J ' u )i< flM <d G R dxd is a strictly positive definite, symmetric matrix and the coefficients 
g K , ft G IN, are random functions depending on t and x such that the mappings Q x [0, T] E5 
(cj,t) i-» g K (u,t, •) are predictable processes with values in certain function spaces. For 
details see Section 12.31 

Equation ([1]) is understood in a weak or distributional sense, i.e. if is a solution of ([1]), 
if for all smooth and compactly supported test functions ip G Cq°(0) the equality 

(u(t,-),p) = (u ,p)+J2 / (a^u^is, -),p}ds + J2 / (^(s,-)^>^< 

holds for all t G [0, T] P-almost surely. Here and throughout the paper we write {u, ip) for 
the application of a distribution u G T>'{0) to a test function <p G Cq°(0). The existence 
and uniqueness of solutions of equation ((T]), respectively equation fl3]) below, within certain 
classes ^ (C?, T) of stochastic processes has been shown in [30]; see also the earlier papers 
by Krylov, Lototsky and Kim, e.g. [29], [31], [32], [35]. Roughly speaking, the classes 
Sjp (O,T) are L p -spaces of functions on Q x [0,T] with values in weighted Sobolev spaces 
Hp d _ p (0) that can be regarded as generalizations of the classical Sobolev spaces with zero 
Dirichlet boundary condition. Again we refer to Section 12.31 for precise definitions. Let us 
remark that in Examples H~T | [TBI and [T9" | illustrating our Besov regularity result in Section [31 
the solution of equation (jlj in the class $)], g (0, T) coincides with the unique weak solution 
with zero Dirichlet boundary condition in the sense of Da Prato, Zabczyk [18], and 
hence can be represented by the well known stochastic variation-of-constants formula 

u(t, ■ ) = e tA u + [ e ( *- s)A G(s) dW s , te[0,T). (2) 
Jo 

Here (e tA ) t >o is the semigroup of contractions on L 2 {0) generated by the partial differential 
operator A = u=l a^ u 9x d2 gx with zero Dirichlet boundary condition considered as an 
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unbounded operator on L 2 (0), (G(i)) t e[o,T\ is an operator-valued process and (Wt)te[o,T\ is 
a cylindrical Wiener process on ^(IN), see Remarks IT31 and [TH in Section [2751 

As already mentioned, our motivation to study the Besov regularity of SPDEs is the 
theme of nonlinear approximation of the solution processes. For deterministic settings, a 
detailed overview of nonlinear approximation and an exposition of the characterization 
of its efficiency in terms of the Besov smoothness of the target functions can be found 
in DeVore [21], see also Cohen [8, Chapters 3 and 4]. Let us consider an example of 
approximation by wavelets in L p (0), the L p -space of real- valued functions on O, p G 
(l,oo). To this end, let {ip\ : A G V} be a wavelet basis on O and let / G L p (0) be 
a target function which we want to approximate by functions /jv G L p (0) belonging to 
certain approximation spaces Sn, where A" is the number of parameters used to describe the 
elements of SV We specify the index set of the wavelet basis by writing V = |J, > , _ 1 V,-; 
the wavelets ipx, A G Vj, j > jo, are those at scale levels j > j respectively, and 
ipx, A G Vj _i, are the scaling functions at the coarsest level jo £ In the case of uniform 
wavelet approximation up to a highest scale level j — 1 + n, n G IN, the approximation 
spaces are 

{3 0-l+ w "| 
5^ Cx ^ x : Cx G R > ^ e V i' i e Oo - 1, • • • , jo - 1 + > , 
j=j -i Aev^- J 

where A" = N(n) = | Uj^ -i ^ ^ * s ^ e car dinality of the set of all indices up to scale 
level j — 1 + n. Let ej^(f) = mif N€SN \\f — /at||l p (o) be the corresponding approximation 
error measured in L p (0). It is well known that — under certain technical assumptions on 
the wavelet basis — the decay rate of ejv(/) is linked to the L p -Sobolev smoothness of the 
target function. More precisely, there exists an upper bound r G IN depending on the 
wavelet basis such that, for all s G [0, r], 

/ e W;(0) e N (f) < C ■ N- s '\ N = N(n), n G IN, 

for some constant C > which does not depend on N. The fractional order Sobolev spaces 
Wp{0) are defined in the next section. One can also show the converse 

3C>0^nem : e N {f)<C-N' s/d , N = N{n) / G W p (O), s' < s. 

If we consider instead best A^-term approximation as a form of nonlinear approximation, 
the approximation spaces are 

Z N = I Y, c ^ : A C V, |A| < iV, c A g E, A G A I, 
I AeA J 

A^ G IN, and in this case the decay rate of the error a^(f) := mff N€ z N \\f — /jv|U„(o) 
is governed by the smoothness of / measured in certain L r (0) -norms, r < p, which are 
weaker than the L p ((9)-norm: For all a G [0,r], 

feB« T (G), - = ^ + - =^ (r N (f)<C-N^ d , New, 
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Figure 1: Linear vs. nonlinear approximation 

illustrated in a DeVore-Triebel diagram. 



B" T (0) being a Besov space as defined in Section Therefore, if the target function / 
belongs to T *(0), 1/r* = a/d + l/p, for some a G [0, r], and if in addition /3 := sup{s G 
R : / G W^(C?)} < a, then the convergence rate of uniform wavelet approximations is 
inferior to the convergence rate of the best iV-term wavelet approximation. The latter can 
be considered as a benchmark for the convergence rate of adaptive numerical algorithms, see 
[2], [12]) [12] • This situation is illustrated in Figure [TJ where each point (1/r, s) represents 
the smoothness spaces of functions with "s derivatives in L T (0)" . Note that the nonlinear 
approximation line {(1/r, s) G [0, oo) 2 : 1/r = s/d + 1/p} is also the Sobolev embedding 
line. For bounded domains, all spaces left to this line as well as the spaces B^. T (0) on the 
line are continuously embedded in L p (0). 

Let us return to equation (jTJ) and assume that the solution u = u(uj,t,x), (u>,t,x) G 
Q x [0, T] x O, vanishes on the boundary dO, satisfying a zero Dirichlet boundary condition. 
It is clear that the smoothness of x h- > u(u,t,x) depends on the smoothness of the mappings 
x i — y g K (u,t,x), k G IN. However, even if the spatial smoothness of the g K is high, the 
Sobolev smoothness of x i— > u(u,t,x) can be additionally limited by singularities of the 
spatial derivatives of u at the boundary of O, due to the zero Dirichlet boundary condition 
and the shape of the domain. Such corner singularities are typical examples for the fact 
that the spatial L p -Sobolev regularity of u may be exceeded by the regularity in the scale 
of Besov spaces B" T (0), 1/r = a/d + l/p. In this paper, we present a result on the spatial 
Besov regularity of the solution u to equation ([T]) which has the following structure: If 

u G L P {Q x [0,T],V, P ® A; W;(0)) 
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and if the functions g K , k G IN, are sufficiently regular, then 

u G L T (Q x [Q,T],V,F®\; B« T (0)) 

for certain a > s and 1/r = a/d + 1/p. Here V is the predictable cr-algebra w.r.t. the 
filtration (J r t)t e [o,T] an d A denotes Lebesgue measure on [0, T]. This result is important for 
the theoretical foundation of adaptive numerical methods for the approximation of u. The 
proof is based on a wavelet expansion of an extension of O 3 x t— > u(u,t,x) to It , which 
allows us to estimate the 5" T (C)-norm in terms of the wavelet coefficients. We apply a 
strategy similar to the one used in Dahlke, DeVore [T3], where the Besov regularity of 
(deterministic) elliptic equations on Lipschitz domains is investigated with the help of an 
estimate of weighted Sobolev norms of harmonic functions. Our substitute for the latter 
is an estimate of weighted Sobolev norms of the solution of equation (JTJ provided by Kim 

There exists an extensive literature on the Besov regularity of SPDEs. In general, 
however, the assumptions on the domain and the scale of parameters considered do not 
fit into our setting. To mention an example, the semigroup approach to SPDEs of Da 
Prato, Zabczyk [18], which is placed in a Hilbert space framework, has been generalized 
to M-type 2 Banach spaces by Brzezniak |3], [5], for the purpose of gaining better Holder 
regularity results. Roughly speaking, the operator A appearing in equation (j2J) is considered 
as the generator of a semigroup on L p (0) for some p > 2, and the stochastic integral in 
([2]) is considered as an stochastic integral in an interpolation space X between L p (0) and 
D(A) C L p {0), the domain of A, realizing a zero Dirichlet boundary condition. If DO is 
sufficiently smooth, then D(A) = W%(0) n W p \0) and X C B s p2 (0) for some s G [0, 2]. In 
this situation, the Sobolev embedding theorem leads to Holder regularity results, and these 
results become better for large p. With the help of a theory of stochastic integration in wider 
classes of Banach spaces, this approach has been generalized in the works of Van Neerven, 
Weis, Dettweiler and Veraar, see, e.g. [20], [37j, [38], [39], compare also Brzezniak, 
Van Neerven [6] . In contrast to these works the problem considered here is of a different 
nature. Firstly, we are explicitly interested in domains with non-smooth boundary. For 

o 

polygonal non-convex domains, it is well known that W%{Q) D W^iO) £ D(A), where 
D{A) := {u e W}{0) : Au G L 2 (0)}, A = A = J^=i & see Grisvard [21], [25], and 
for more general Lipschitz domains see Jerison, Kenig [28J. Secondly, we are interested 
in the special scale B" T (0), 1/r = a/d+ 1/p, r > 0, p fixed, including in particular 
spaces which are no Banach spaces but quasi-Banach spaces. The parameter r decreases 
if a increases and B™ T {0) fails to be a Banach space for r < 1. While our methods work 
in this setting, any direct approach requires (at least!) a fully-fledged theory of stochastic 
integration in quasi-Banach spaces which is not yet available. 

Let us emphasize that our result can be extended to more general linear equations of 
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the type 



fl,V=l K=l \(1=1 

u(0, ■) = u , 



(3) 



including, in particular, the case of multiplicative noise. Here the coefficients a^ u , U\ c, 
cr M/t , rf and the free terms / and g K are random functions depending on t and x. This 
extension is possible because one of our main tools, the weighted Sobolev norm estimate 
of Corollary [12] in Section 12. 3[ holds for equations of type (pQ) as well as for equations of 
type ([3]). Since this mainly adds notational complications, we will focus on equation ([1]) 
and refer to Appendix [B] for a short account of how to treat equations of type (|3j). 

The paper is organized as follows: In Section 2 we collect the notations, definitions 
and preliminary results needed later on. Some general notations are introduced in Section 
12.11 Section 12.21 provides the necessary facts on Besov spaces and wavelet decompositions. 
In Section 12.31 a short introduction to the general L p -theory of SPDEs on Lipschitz do- 
mains due to Kim [30J is given, including definitions of the already mentioned spaces 
Hp g __ p (0), f)p 8 (0,T). Finally, in Section |3] the Besov regularity result (Theorem IT5|) is 
stated and proved, and some concrete examples for an application of the result are given. 



2 Preliminaries 

2.1 Some notations and conventions 

In this and the next subsection O C R d can be an arbitrary (not necessarily bounded) 
Lipschitz domain. A domain is called Lipschitz if each point on the boundary dO has a 
neighbourhood whose intersection with the boundary — after relabeling and reorienting the 
coordinate axes if necessary — is the graph of a Lipschitz function 

By T>'(0) we denote the space of Schwartz distributions on O. If not explicitly stated 
otherwise, all function spaces or spaces of distributions are meant to be spaces of real- valued 
functions or distributions. If / G V(0) is a generalized function and a = (ckj, . . . , aj) G 
INq is a multi-index, we write D a f = a f ^ dx <* d f° r the corresponding derivative w.r.t. 
x = (xx, . . . , Xd) G O, where \a\ — ai + . . . + As in equations ([T]) and ()3]) we also use 
the notation f X)lXu = f Xfi = |£. For m G M„, D m f = {D a f : \a\ = m} is the set 

of all m-th order derivatives of / which is identified with an Rv m /-valued distribution. 
Given p G [1, oo) and m G INq, W™(0) denotes the classical Sobolev space consisting of all 
(equivalence classes of) measurable functions / : O — > R such that ||/||vy m (e>) = ||/||l p (0) + 
\f\w P (o) = (Jo Ifix^dx) 1 /? + E\ a]=m Uo \D a f{x)\vdxy/v is finite. Fo'r p G (l,oo) and 
s G (m,m + 1), m G 1N , we define the fractional order Sobolev space Wp(0) to be 
the Besov space Bt, JO) introduced in the next subsection. (This scale of fractional order 
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Sobolev spaces can also be obtained by real interpolation of Wp(0), n G 1N . One can 
show that W£(0) = S£ 2 (C>) for all n G IN and W£{0) C B^ p (0) for all n G IN, p > 2, 
see, e.g. Triebel [HJ Remark 2.3.3/4 and Theorem 4.6.1.(b)] together with Dispa [23].) 
Given any countable index set J, the space of p-summable sequences indexed by J is 
denoted by £ p = £ P (J) and | ■ \g p is the respective norm. Usually we have £ p = £ p (ti$) but, 
for instance we may also use the notation \D m f(x)\ p = J2\ a \=m \D a f(x)\ p for / G W™{G). 

Given a distribution / G D'(0) and a smooth and compactly supported test function 
ip G Cq°(0), we write (/, ip) for the application of / to <p. If if is a Hilbert space, then 
( • , • )h denotes the inner product in H. Given another Hilbert space U, we denote by 
L(Hs)(if, U) and L( nuc )(ii", U) the spaces of Hilbert-Schmidt operators and nuclear operators 
from H to U respectively, see, e.g. Pietsch Sections 6 and 15] or Da Prato, 
Zabczyk [TSJ Appendix C] for definitions. We also abbreviate L^(H) = L^(H,H) 
and L( nuc )(i?) = L( nuc )(i7, H). J\4^°(H, (J r t)) is the space of continuous, square integrable, 
if- valued martingales with respect to the filtration (J-t)te[o,T\- For Q x [0,T] we use the 
shorthand notation an d 

V = a({]s,t] x F s : < s < t < T, F s G F a } U {{0} x F : F G ^o}) 

is the predictable cr-algebra. P <8> A is the product measure of the probability measure P 
on (Q, J 7 ) and Lebesgue measure A on ([0, T], B([0, T])), where £>([0,T]) denotes the Borel 
a-algebra on [0, T]. Given any measure space (A, A, m), any (quasi-)normed space B with 
(quasi-)norm \\-\\b and any summability index p > 0, we denote by L P (A, A, m; B) the L p - 
space of all strongly measurable functions u : A B whose (quasi-)norm ||m||l p (a,A"i;B) : = 
(J A \\u(z)\\ p B m(dz)) l ^ p is finite. 

All equalities of random variables or random (generalized) functions appearing in this 
paper are meant to be P-almost sure equalities. Throughout the paper, C denotes a positive 
constant which may change its value from line to line. 



2.2 Besov spaces and wavelet decompositions 

In this section we give the definition of Besov spaces and describe their characterization in 
terms of wavelets. Our standard reference in this context is the monograph of Cohen [S]. 
For a function / : O — > R and a natural number n G IN let 



n n / \ 

i=0 j=0 ^ ' 



(-l)"-'f(x + jh) 



be the n-th difference of / with step h G R d . For p G (0, oo) the modulus of smoothness is 
given by 

u n (tj) p := sup \\A n h f\\ Lp(oh t>0. 

\h\<t 

One approach to introduce Besov spaces is the following. 
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Definition 1. Let s,p,q G (0, oo) and n G IN with n > s. Then B^ q {0) is the collection 
of all functions / G L p (0) such that 



|/|i»|,,(0) ■= [t S U n (tJ) p \ j | < X. 

These classes are equipped with a (quasi-)norm by taking 

II / ||b« 9 (0) : = II / lUp(O) + I / ls|,,(o) • 

Remark 2. For a more general definition of Besov spaces, including the cases where 
p, q = oo and s < see, e.g. Triebel [35] . 



We want to describe -Bp g (R d ) by means of wavelet expansions. To this end let p be a 
scaling function of tensor product type on R d and let ipi, i — 1, . . . , 2 d — 1, be corresponding 
multivariate mother wavelets, such that, for a given r G IN and some N > 0, the following 
locality, smoothness and vanishing moment conditions hold. For all i = 1, . . . , 2 — 1, 

supp (p, supp tpi C [-JV, iV] d , (4) 
<p, ^ G C r (R d ), (5) 

J x a i)i{x) dx = for all a G JN„ with |a| < r. (6) 



We assume that 



{Vk,^,k : (*,J,fc)G{l,---,2 d -l}xlNoxZ d } 



is a Riesz basis of L2(R d ), where we use the standard abbreviations for dyadic shifts and 
dilations of the scaling function and the corresponding wavelets 

ip k (x) := <p(x - k), x G R d , for k G Z d , and (7) 

ipi,j,k(x) := 2 jd/2 ipi(2 j x - k), x G R d , for (z, j, jfe) G {1, • • • , 2 d - 1} x IN Q x Z d . (8) 

Further, we assume that there exists a dual Riesz basis satisfying the same requirements. 
More precisely, there exist functions tp and i — 1, . . . , 2 d — 1, such that conditions (jlj), 
(J5J and (jSJ) hold if <p and ^ are replaced by <p and ipi, and such that the biorthogonality 
relations 

(<Pk, i>i,j,k) = (^i,j,k, <Pk) = , (tpk, fi) = 8k,i, {4>i,j,k, i>u,v,e) = $i,u 8j,v h,£ , 

are fulfilled. Here we use analoguous abbreviations to (j7j) and (jSJ) for the dyadic shifts 
and dilations of <p and ipi , and 5^ denotes the Kronecker symbol. We refer to Cohen [SJ 
Chapter 2] for the construction of biorthogonal wavelet bases, see also Daubechies 
and Cohen, Daubechies, Feauveau [IT] . To keep notation simple, we will write 

Aj,k, P ■= 2 jd{1/p - 1/2) A,j,k and := T?^- 1 ' 2 ^^ 
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for the Lp-normalized wavelets and the correspondingly modified duals, with p' := pj (p — 1) 
if p G (0, oo), p 7^ 1, and p' := oo, 1/p' := if p = 1. 

The following theorem shows how Besov spaces can be described by decay properties 
of the wavelet coefficients, if the parameters fulfil certain conditions. 

Theorem 3. Let p,q G (0, oo) and s > max{0, d (1/p — 1)}. Choose r G IN such that 
r > s and construct a biorthogonal wavelet Riesz basis as described above. Then a locally 
integrable function f : R d — >■ R is in the Besov space Bp q (R d ) if, and only if 

2 d -l 

f = E ^' ^ ^ fc + E E E ^' ^> fc >p') ^i-k,P ( 9 ) 
(convergence in T>'{K d )) with 

2 d — 1 

( E k/^>i p ) 1/p + (E E 23sq { E \(f^M)\ p ) q y q < °°> 

and ffTOj) zs an equivalent (quasi-)norm for Bp q (R d ). 

Remark 4. A proof of this theorem for the case p > 1 can be found in Meyer [361 §10 
of Chapter 6]. For the general case see for example Kyriazis [34] or Cohen [8, Theorem 
3.7.7]. Of course, if ffTOj) holds then the infinite sum in converges also in Bp q (R d ). If 
s > max {0, d (1/p — 1)} we have the embedding B^ q (R d ) C L u (R d ) for some u > 1, see, 
e.g. Cohen [SJ Corollary 3.7.1]. 

Let us now fix a value p G (1, oo) and consider the scale of Besov spaces S^ T (R d ), 
1/r = s/d + 1/p, s > 0. A simple computation gives the following result. 

Corollary 5. Let p G (1, oo), s > and r G R stzcn taat 1/r = s/d + 1/p. Choose r G IN 
siic/i t/iat r > s and construct a biorthogonal wavelet Riesz basis as described above. Then 
a locally integrable function f : R d — > R is in the Besov space B^ T (R d ) if, and only if 

2 d -l 

/ = E $k)vk + E E E ^^p') ^jap ( n ) 

fcgz d «=i iew fcez d 

(convergence in T>'(Jl d )) u>rfn 

2 d — 1 

( E i(/.^)r) 1/T +(EE E \(f^M)\ T ) lT < (12) 

feez d i=l jelNo k&A 

and (|T2|) an equivalent (quasi-)norm for B* T (R d ). 
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2.3 SPDEs on Lipschitz domains and weighted Sobolev spaces 

From now on, let O C R d be a bounded Lipschitz domain and d > 2. 

We have already mentioned corner singularities as typical examples where the regularity 
of a function on O C R d in the Besov scale B" T (0), 1/r = a/d+l/p, a > 0, can exceed the 
regularity in the Sobolev scale 1^(0), s > 0. This reflects the sparsity of the large wavelet 
coefficients of such a function (given a wavelet basis on the domain O). A general way to 
deal with smoothness regardless of certain singularities at the boundary is to use weighted 
Sobolev spaces, where the weight function is a power of the distance to the boundary. 
The Lp-theory of SPDEs on Lipschitz domains by Kim [30J is based on spaces of this 
type, namely the weighted Sobolev spaces H^ g (0), p G (l,oo), 9, 7 G R, introduced in 
Lototsky |35J. They are defined in terms of the Bessel-potential spaces 

H~<(R d ) = {ue S'(R d ) : \\u\\ H , (Rd) = ||(1 - Ay/ 2 u\\ Lp(Rd) < 00}. 

Here, <S'(R d ) C V(R d ) is the space of (real valued) tempered distributions and (1 — A) 7 / 2 : 
S'(R d ) S'(R d ) is the pseudo-differential operator with symbol R d 3 £ h-> (1 + |£| 2 ) 7/2 , 
i.e. (1 — A) 7//2 m = J-" _1 ((l + l^l 2 ) 7 / 2 ^?/), where J 7 denotes the Fourier transform on the 
(complex valued) tempered distributions. 

For iG0we write p(x) := dist(x,dO) for the distance between x and the boundary 
of the domain O. Fix c > 1, k$ > and for n G Z consider the subsets O n of O given by 

O n := {x G O : c- n - fco < p(x) < cT n+fe °}. 

Let Cn, n G Z, be non-negative functions satisfying £ n G C£°(O n ), ^ ng ^Cn(a ; ) = 1 and 
|-D m Cn(^)| < C ' c mn for all n G Z, m G JN , x G 0, and a constant C > that does not 
depend on n, m and x. The functions ( n can be constructed by mollifying the indicator 
functions of the sets O n , see, e.g. Hormander [271 Section 1.4]. If O n is empty we set 
( n = 0. For u G V{0) ( n ii is a distribution on with compact support which can be 
extended by zero to R d . This extension is a tempered distribution, i.e. ( n u G iS'(R d ). 

Definition 6. Let ( n , n G Z, be as above and p G (1, 00), 9, 7 G R. Then 

According to Lototsky [35J this definition is independent of the specific choice of 
c, ko and £ n , n G JN , in the sense that one gets equivalent norms. If 7 = m G INq then the 
spaces can be characterized as 

H%{0) = L p , e (0) := L p (0, P (x) e ~ d dx), 

H™ e (0) = {u : p H D a u G L pfi iO) for all a G < with |a| < m} , 
and one has the norm equivalence 

< / |^) |a| ^ a w(^)r^) <, ~ d da; < C||n||^ (0) . (13) 
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Analogous notations are used for t 2 = ^2(lN)-valued functions g = (g K ) K e$!- For p G 
'1, oo), 9, 7 G R and ( n , n G Z, as above 

iy;(R d ;£ 2 ) := G (5'(R d )) ]N : (1 - A) 7/ V e ^ P (^ d ) for all k G IN and 

IMk(W 2 ) == || I ((1 - A ) 7/ V) KeW |, 2 || Lp(Rd) < oo}, 



Remark 7. (a) One can consider the spaces H^ e {0) as generalizations of the classical 
Sobolev spaces on O with zero boundary conditions. For 7 = m G INo we have the identity 

and the norms in both spaces are equivalent, see Theorem 9.7. in Kufner [33]. Here 
W™{0) is the closure of C °°(£>) in the classical Sobolev space W™{0). 

(b) Note that, in contrast to the spaces W£((D) = Bp JO), s G (m,m + 1), m G 
No, which can be regarded as real interpolation spaces of the classical Sobolev spaces 
W™{0), m G IN (see, e.g. Triebel 03 Section 1.11.8] and Dispa[23J), the spaces 
Hp 6 (0), 7 G (m,m + 1), m G IN , are complex interpolants of the respective integer 
smoothness spaces (cf. Lototsky [35] Proposition 2.4]). 

We can now define spaces of stochastic processes and random functions in terms of the 
weighted Sobolev spaces introduced above. 

Definition 8. For 7, 9 G R and p G (1, 00) we set 

T) := L p (Or, V, P ® A; , 
Bj, (0, T; 40 := L p (fi T , P <g> A; fl^ (0; £ 2 )) , 

^(O) :=L„(n,.Fo,P; ^lp(O)), 

and for p G [2, 00), 

{00 
u G H^_ p (0, T) : «(0, • ) G E£,(0) and <fc = / dt + £ <f <K 

for some / G H^(0,T), £ G EjftO, T;4)}, 

equipped with the norm 

\Wh; i9 (o,T) ■■= \\u\\ml e _ p (o,T) + II/IIh^ p (o,t) + IMIh^co.t^) + H M ( ' -)lli£ 9 (o)- 
The equality du = f dt + Yl™=i g K dw£ above is a shorthand for 

(u(t, .),<p)= (u(0, -),<p)+ (/( s > -),<P)d* + Y, ^ s > -),<p) dw: (14) 

^0 K=1 Jo 

for all y> G Cg°(<9), t G [0,T]. 
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Remark 9. (a) If p G [2, oo), then the sum of stochastic integrals in (1141) converges in the 
space A4y C (R, (J 7 *)) °f continuous, square integrable, R-valued martingales w.r.t (J~t)te[o,T\- 
For the convenience of the readers we include a proof in Appendix |A] 

(b) Using the arguments of Krylov in [31] Remark 3.3], we get the uniqueness (up to 
indistinguishability) of the pair (f,g) G E£? +p (0,T) x Hj/(0,T;4) which fulfils (JH}. 
Consequently, the norm in $)p d (0,T) is well defined. 

Definition 10. We call a predictable V'(0)-vaAued stochastic process u = (u(t, -))*e[o,T] 
a solution of equation (JT]) if it is a solution of equation f JT4|) where / is replaced by 
Y?p,v=i a ^ u x»x v and w(0, ■ ) = u . 

The next result is taken from Kim [3T?] . 

Theorem 11. Lei p G [2, oo) and 7 G R. There exists a constant k G (0,1), depend- 
ing only on d, p, {a^ v )x<^v<d and O, such that for any 6 G (d — k , d — 2 + p + k ), 
g G Hp ~^{0, T; £2) and Uq G Up 9 (0) the equation (JTJ) /ias a unique solution in the class 

For this solution 



where the constant C depends only on d, p, 7, 6, {o> IMV )i<n ) i/<di T and O. 

We will need the following straightforward consequence of this Theorem [TTJ Recall that 
if m G IN and / G V(0) is sufficient regular, then \D m f\ ip stands for (£ H=m \D a f\ p ) 1/p , 
the (pointwise) £ p -novm of the vector of the m-th order derivatives of /. 

Corollary 12. In the situation of Theorem [771 with 7 = m G IN, the following inequality 
holds for every r G [0,p\. 

\p m ~ s \D m u(u,t, •)kllMO)* p ( dw ) ^ C (II^IIh™- 1 ^,^) + NI|t/™(0)) , 

where 5 = 1 + — . 

p 

Proof. Theorem [11] implies, in particular, that 

IHIh™ 9 _ p (c>,t) < C (\\g\\ M ™~i( ,T) + \\ u o\\u™ e (o) 

and we have 

r-T 




Q JO 



Mle™ (ot)= / OlljU, (OT) diP(dw) 



fc=0 
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n Jo 



>C / \\p m -°\D m u{u,t,.)\ t J L , 0) dtI>{du) 
Jn Jo 

with 5 = 1 + G ((2 — K )/p , (p + Kq)/p). Now let r G [0,p\. Jensen's inequality for 
concave functions, see, e.g. Schilling [121 Theorem 12.14], yields 

T \\p m - s \D m u(u,t, .)\ e Jl p(0) dtF(duj) < C(T) (lbll^- 1(0jT) + ll«o||^ (0) N ^ 

< c (lloWm^- 1 (o,T) + \\ u o\\u™ e (o)^ ■ 

In the last step we have used the fact that all norms on R 2 are equivalent. □ 

Remark 13. Consider the Hilbert space case p = 2 and assume g G EL] g (0, T; I2). The 
expression Y^=\ Jo 9 K ( S > ' ) can ^ e considered as an H] e (C)-valued stochastic integral 
f G(s) dW s with respect to a cylindrical Wiener process (Wt)te[o,T\ 011 ^2 whose coordi- 
nate processes are (w£)te[o,T\, ^ G IN. (See, e.g. Da Prato, Zabczyk [18J or Peszat, 
Zabczyk jlQ] for stochastic integration w.r.t. cylindrical processes.) Here (G(t))te[o,T] is a 
stochastic process in the space of Hilbert-Schmidt operators L(hs)(-^2, H^ei®)) defined by 

G(u, t):£ 2 ^ Hl $ (0), (x K ) KeW 1 y 9*( u > £ > • ^ £ ) G 

and it is an element of the space L 2 (^t! Lms) (^2, H^^O))). Indeed, for fixed (u,t) G Qt 
we have 

\\G(^m% HS)i e 2 ,Hi em = E H^( w '*» Oll^o) 

= E E Cne !K-n(c" • )<?>, f, C" • )||^ (Rd) 

2 



Ec n 1C-n( C "-)^(^^c" 



l^(R d ;£ 2 ) 



by Tonelli's theorem, so that 

||G|| L2 (n T ;i (HS) (fe;ifT e ( 0))) = |b||Hj 9 (0,T;&)- 

As a consequence, equation can be rewritten in the form 

d 

du= E 0""^^ * + dAf t , u(0, • ) = u , (16) 

where (M t ) te [ 0j r] G M 2 /(H^ e (0), (Ji)) is the # 2 7 e (0)-valued, square-integrable martingale 
given by 



t ■- 



/ G(s)dW s , t G [0,T]. 
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Remark 14. In Examples IT?] ITBl and [T9l below the solution u of equation (P) in ^ e {0,T) 
as given by Theorem [TT1 coincides with the weak solution of equation ( fl6|) with zero Dirichlet 
boundary condition in the sense of Da Prato, Zabczyk [T8"] . 

In the examples we consider equation (FTO]) driven by certain Wiener processes (M t ) te [ 0)T ] 
in L 2 (0) with u G (7 2 2 2(0), d = 2 and the solution u is in the class ^22(^5^) c 
1120(0, T). (Strictly speaking, in Example IT9l (MA^ht] is not a Wiener process, but it is 
one conditioned on the family of random variables Y\, A G V.) Thus, by Remark [7] (a) we 

o 

know that u is an element of L 2 (VLt] W$ \(0)). Let us now introduce the operator 

(A, D(A)) := ( a " V o^T . {« e #2(0) : E e L 2 (C) 

and consider the equation 

du(t, ■) =Au(t, ■) dt + dMt, u{0, ■) = u eL 2 (O), te[0,T}. (17) 

A weak solution of equation ffTTj) in the sense of Da Prato, Zabczyk [18J is an L 2 (0)- 
valued predictable process u = (u(t, ■ ))te[o,T\ with P-almost surely Bochner integrable 
trajectories t !->■ u(u,t, ■ ) satisfying 

■ ), C)£2(0) = ( u o, C)l 2 (o) + / («(s, ■ ), A*C)l 2 (0) ds + (M t , C)i a (o) (18) 

./o 

for all t G [0,T] and £ G .0(^4*). It is given by the variation of constants formula 

u (t, ■ ) = e tA u + [ e {t - s)A dM s , t G [0, T], 
Jo 

where (e M ) t > is the contraction semigroup on L 2 {0) generated by A. 

It is clear that the solution u G ,£j| 2 ((D, T) given by Theorem UTI satisfies 

(u(t, ■ ), <p) L2{0) = (u , <p) L2io) + / (A^uis, ■ ), A 1 ' 2 ?)^) ds + (M t , V ) L2(0) (19) 

Jo 

for all t G [0,T] and if G Cq°(0). Note that the operator A is self-adjoint because the 

c 

coefficients a^, 1 < /i, f < d, are constants. Since every ( G D(A*) = D(A) C W 7 ^ f(C) is 
the limit in W 2 1 (0) of a sequence of test functions ((pk)ke¥i C Cq°(0), one can go to the 
limit k — > 00 for ip = ifk in (|T9|) to obtain equation ffTSl . 
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3 Besov regularity for SPDEs 

In this section we state and prove our main result. We give some concrete examples to 
illustrate its applicability. The result is formulated in terms of the L T -spaces 

L T (n T ; B S T>T (0)) = L T (n T , V, P ® A; fl' T (0)), t G (0, oo), s G (0, oo), 

and the spaces introduced in the last section. 

Theorem 15. Let p G [2, oo) and g G H^O, T; £ 2 ), u G U^ e (0) for some 7 G IN and 
6 G (d — Kq, d — 2 +p + Kq) with k = K (d,p, {a^), O) G (0, 1) as in Theorem [771 Let u be 
the unique solution in the class Sj 1 {O,T) of equation (JTJ) and assume furthermore that 

d-9\i 



u G L P (Q T ; B s p>p (0)) for some s G (o, 7 A (l + — — ) . (20) 
Then, we have 

I cm 1 / id s 

ueL T (Q T -B: T (0)), - = -7 + - for all a6 0, 7 A — 

ran V a — 1, 



P 

and i/ie following estimate holds 

\\u\\ Lt (q t ;B ?iT (o)) < C (iMI IT" 1 (0,T;4) + H M o||cryo) + ||w||i s (n T ifli, p (0))) • ( 21 ) 
i/ere ine constant C depends only on d, p, 7, a, s, 0, {d^ v )i<n,u<d-, T and O. 



Remark 16. Since the constant k = K Q (d,p, (a^), O) is greater than zero, we can always 
choose 9 = d. In this case, we know from Theorem [TT] that for each 7 G IN we have a unique 
solution u in the class Sjl d (0, T), provided the free term g and the initial condition uq are 
sufficiently regular. In particular, we get 

u G Ml d _ p (0, T) = L P (Q T , V, P ® A; H^(0)) C L p (fi T ; Wj(O)) C L p (fi T ; 5^(0)). 

Thus, the additional requirement ( 120]) is fulfilled with s = 1. Since O is an arbitrary 
bounded Lipschitz domain it is in general not clear if u belongs to L p (Qt', W£(0)) for all 
s < 2, compare Example [T7I below. 

However, if 7 > 2 our result shows that we obtain higher regularity than s = 1 in the 
nonlinear approximation scale, namely 

u G L T (n T ; B° T (0)), - = 3 + -, for all a < 

t d p d — 1 
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Proof of Theorem [T5[ We fix a and r as stated in the theorem and choose a wavelet Riesz- 
basis 

{<p k , V M > = (i,3, k) e {1, • ■ • , 2 d - 1} x ]N x Z d } 

of L2(R d ) which fulfils the assumptions from Section I2~2l with r > 7. Given (j, k) G INq x Z d 
let 

g jifc :=2-^ + 2^[-iV,iV] d , 

such that supp ipij )k C Qj >k for alH G {1, . . . , 2 d — 1} and supp C Qo,fc f° r a U A; G Z d . 
Remember that the supports of the corresponding dual basis fulfil the same requirements. 
For our purpose the set of all indices associated with that wavelets that may have common 
support with the domain O will play an important role and we denote them by 

A := {(hj, k) G {1, . . . , 2 d - 1} x 1N x Z d | Q jjk nO^}. 

In particular, we will also use the following notation: 

r := {k G Z d : Q ,fc H O £ 0}. 

Due to the assumption u G L p {VLt] B^ p {0)) we have u(u,t, •) G B* (O) for P <g) A- 
almost every (u, t) G fiy. As O is a Lipschitz domain there exists a linear and bounded 
extension operator £ : Bt, JO) — >■ 5p p (R d ), i.e. there exists a constant C > such that for 
P <g> A- almost every (co,t) G O r : 

Su(w,t, -)\ = u(u,t, •) and \\Su(u,t, ■ )\\ B ^ p (M,d) < C\\u(u,t, -)b|, p (o)> 

see, e.g. Rychkov [¥T] . In the sequel we will omit the £ in our notation and write u instead 
of £ u. 

Theorem [3] tells us that for almost all (00, t) G Qt the following equality holds on the 
domain O 

uiu, t, ■ ) = ^2(u(u, t, ■ ), (p k )<f k + ^2 *> ' )' $i,j,k,p')i>i,j,k, P , 

fcer (i,j,k)eA 

where the sums converge unconditionally in Sp p (R ). Furthermore, cf. Corollary [5j we get 
for P ® A-almost all (to, t) G Vt T 

\\u(u,t, -Wb^o) <c(j^\{u{uj,t, ■),^ k )\ T + K"M, •),^,i,fcy)r)- (22) 

fcer («,i,fc)6A 

Hence, it is enough to prove that 

/ / £|<u(a;,t, -j-^r^PW < CI^IIm^bi ,(©)) (23) 
Jn Jo Ul -r 
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and 



I I J2 \(u(u,t,-),ip i j Ap ,)\ T dtF(du J ) 
JnJo (i,j,k)eA 

< c (IIpIIbj-^©.^) + IKHt/yo) + \\u\\L p (n T ;B^ p (o))^) 



(24) 



We start with (T231 . The index set T introduced above is finite because of the bound- 
edness of the domain O, so that we can use Jensen's inequality to get for P ® A-almost all 
(co,t) G Q T 

1 t, .), ip k ) r < c ( ( J2 1 M w > *» ■). i p ) 1/P ) < <? n«(w. *, o ii^(o)- 

In the last step we used Theorem [3] and the boundedness of the extension operator. Inte- 
gration with respect to P <8> A and another application of Jensen's inequality yield (I2"3"j) . 
Now let us focus on the inequality (j24p . To this end, we introduce the following notation 



Pj,k 


:=dist(Q j>fc ,<9(9) 


= inf p(x), 




A, 


:= {(i,l,k) G A : 






:= {{i,j,k) G A, 


: ^2 _i < p j>k < (m + 1)2" 


A? 


■ \ \ A, . 




A 


== U A?, 

j€W 





where j, m G 1N and k G Z d . We split the expression on the left hand side of (I24p into 

/ / K u ( w '*> ' ),4>i,j,k, P >)\ T dtF(du) 

Jn Jo ., v ,„ 




a Jo 



\(n(u,t, ■),^i, j ,k, P '}\ T dtF(du) =: I + II (25) 

(i,j,k)eA\A° 



and estimate each term separately 

Let us begin with I. Fix k) G A and (co,t) G Qt such that 



o 



|p(x) 7 s \D' y u(u,t,x)\i p \ p dx < 



oo. 



By Corollary [15] this holds for P ® A-almost all (tu, £) G fir- By a Whitney-type inequality, 
also known as the Deny-Lions lemma, see, e.g. DeVore, Sharpley [221 Theorem 3.4], 
there exists a polynomial Pj^ of total degree less than 7 such that 



\u 



(U,t, •) - Pj,k\\L p (Q j}k ) < C2 n \u(u,t, • 



17 



where the last norm is finite since pj^ = dist(Qj,fc, dO) > 0. Since ?pij,k,p' is orthogonal to 
every polynomial of total degree less than 7, one gets 

I (u(u, t, ■ ),$i,j,k, P >)\ = I t, •)- Pj,k, i>i,j,k, P ') I 

< \\u(u,t, ■ ) - Pj,k\\Lp(Q jlh ) \Wi,jAp'\\L p ,{Q jlk ) 



<C2-^\u(u,t,.)\ w ^ h) 



i/p 



< ci '-,> jt ;{ / \ P ( x y- s \DMu,t,x)\ ip \ p dx) 
■■C2 r n jJ : l > J jU.i). 

Fix j G INo- Summing over all indices k) G A° and applying Holder's inequality with 
exponents - > 1 and one finds 

r T p—T 

E \Hu,t, -),a m )\ t <c e 2-">5;rV J , t (w,«) r 

(ij,fc)eA° (ij,k)6A5 

<c( £ /*,»M)*)*( E 2^^)^.(26) 

Since any x G (9 lies outside of all but at most a constant number C > of the cubes Q/,fc, 
A; G Z d , we get the following bound for the first factor on the right hand side 

T_ T 

( E H k ^^) p \ = { e / i^) 7-s p 7 «(^*^)kr^) p 



<C\\p^ s \DMu,t,-)U P \\ T L p{ o)- ( 27 ) 

In order to estimate the second factor in (I2"6T) we use the Lipschitz character of the domain 
O which implies that 

I Aj, m | < C2 j ^ for all j, m G IN . (28) 

The constant C > does not depend on j or m. Moreover, the boundedness of O yields 
Aj, m = for all j, m G !N"o with m > C2? . Consequently, 

p — T (^03 p—T 



(29) 



2 Pi,r ) <(E E 2p ~ T Pd 

(i,i,fc)eAP y V "i=l (tJ,fc)eA il m 



V »V1 — 1 / 



<C(2^- 1 -F^) + 2 J (' i -^) 



p — T 
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Now, let us sum over all j G 1N and integrate over with respect to P <8> A on both sides 
of the inequality (126]) . By using ( 12"9~j) and (127)1 from above and Corollary [T2] we get 



/ / E K u ( w >*> •)>v , ij,fc,p')r^i > (^) 

(ij,A)eA° 

/ JSl Jo 



One can see that the sums on the right hand side converge if, and only if, a G (0, j A sjzj) • 
Finally, 

// E |<w(^,*, • ),^ij,fcj/>r*P(dw) < C' (llfl'llHj-i^r^) + llwolltg^co)) • (3°) 
' /<! (M,fc)eA° 

Now we estimate the term II in (125|) . First we fix j G INo and use Holder's inequality 
and (J2ED to get 

T 

(i,j,fe)eAj,o (i,j,fc)eAj,o 
Summing over all j G IN and using Holder's inequality again, yields 

E \( u ( u > t >-)>$ij,w)\ T = E [ Yl | <w(w, *, o > ^tj.fcj/) r 

(tj,fc)eA\A° j'eWo (i,j,k)eA j:0 

jGMo (i,j,k)eAj, 

< c (yj 2 ^( I ^-)^) z? (5: yj a^K^,*,.),*^)!')*- 

Using Theorem [3] and the boundedness of the extension operator, one gets for P ® A-almost 
every (w,t) G Or that 

£ \(u(u,t,.U iM )\ T < C\Hco,t,-WB U o)(j2 2^ I ^ li - sp )-) P • 
(ij,k)eA\A° ^jeiNo ' 

The series on the right hand side converges if and only if a G (0, Sjzt)- But this is part of 
our assumptions, so that for P <8> A-almost every (u,t) G Qt 

\{ U {^^)^i M )\ T < C\\u{u,t,.)\\ T B s p{0 y 

(ij,k)eA\A° 
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Let us integrate over Q T with respect to P <8> A and use Jensen's inequality to get 




^2 \(u(w,t, •),Vi,i,fe,p'}| r c?*I > (^) < C / / \\u(u,t, ■) ||fl* p (o) dtF(dw 



n Jo (i,i,fc)eA\A° J>) 



< C 




in Jo 

Because of ([30]) this proves Q2H). Now and Q22|) finish the proof. □ 

Next, we give some examples for an application of Theorem IT5l We are mainly interested 
in the Hilbert space case p = 2. 

Example 17. Let us first consider equation (TTJ) in the form ffT6l) where the driving process 

o o 

(M t ) t€ [o : T] is a Wiener process in W 2 l (0) with covariance operator Q 6 L( nuc )(W 2 (O)). It 
can be represented as a stochastic integral process (/J G(s) dW s ) t ^[o,T] w.r.t. the cylindrical 
Wiener process (W t )te[o,T] on £ 2 by defining the integrand process (G(t))te[o,T] in the space 

o 

of Hilbert-Schmidt operators L(hs)(^2> W% (O)) as the constant deterministic process 



G(u,t) : £ 2 -+ W}(0), (x K ) Km M- V m V^x"e„, (w,t) e fi T , (31) 

o 

where (e K ) K£ ]N is an orthonormal basis of W 2 (C?) consisting of eigenvectors of Q with 
positive eigenvalues (X K ) K ^. 

This corresponds to defining g = (g K ) Ke -^ in equation ([TJ) by 



g^co, t,-):= VA K e K , k G IN, (w, t) G fi T . (32) 
It is easy to see that g is an element of T; £ 2 ). By definition 



= T " E A « E c "1C-n(c" • )e K (c n 011^) 



(33) 



Using the norm equivalence ( lT3"j) . one has 

imi^ow,) < cr 2 E A * E II^'^IIL^) < ct2 E a - E P Qe «iiL(o) 

k£K |a|<l K6K |q|<1 

= CT 2 ^A K < oo. 
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Figure 2: Besov regularity in the scale B^ T (0), 1/r = a/2 + 1/2, vs. Sobolev regularity 
of the solution, illustrated in a DeVore-Triebel diagram. 

Thus, in a 2-dimensional setting, Theorem [TTJ with d = 9 = 7 = 2 tells us that for every 
initial condition u G U 22 {0) — L 2 (Q, J-q, P; -ff 2 2^)) equation (pQ) has a unique solution 
u in the class i}| 2 ((!}, T) c H| o (0, T) = -^(^t; #f o(^))- As a trivial consequence, 

u G L 2 {n T - W 2 X {0)) = L 2 (Q T - B\ 2 (0)) 

because we have the equality 

H 2 2fl (0) = {ue V\0) : p^ 1 D a u G L 2 {0) for all a G INjj with |a| < 2}. 

o 

(In fact, according to Remark [7] we even know that u G L 2 (Q T ; W 2 \0)).) 

Note that in general u does not belong to L 2 (Qt',W 2 (0)) for all s < 2. Since O is 
an arbitrary bounded Lipschitz domain, certain second derivatives might explode near the 
boundary and the norm \\u(co,t, -)llw 2 (0) as well as \\u(u,t, ■ )\\w£(o), where s G (1,2), 
might not be finite. If O is a polygonal domain, one can derive an explicit upper bound 
for the regularity in the Sobolev scale L 2 (fl T ; W 2 (0)), s > 0. Adapting techniques used 
in Grisvard [21], [25] to our stochastic setting, one can show that u ^ L 2 (Q T ; W 2 (0)) if 
s > 1 + 7r/7o, where 70 is the measure of the largest interior angle at a corner of dO. 

However, in the considered situation Theorem [15] with s = 1 states that we have 

ueL T (n T ;B« T (0)) 

for every a < 2 and 1/r = a/2 + 1/2. This constellation is illustrated in Figure EJ where 
each point (1/r, s) represents the smoothness spaces of functions with "s derivatives in 
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L T (0) n . Based on the knowledge that u G L 2 (Q T ; W 2 \0)) and u G L T (Q T ; B% T (0)) for 
all a < 2, 1/r = a/2 + 1/2, interpolation and embedding theorems yield that u also 
belongs to each of the spaces L T {VL T \ B* (O)), < r < 2, s < (1/2 + 1/r) A 2. This is 
indicated by the shaded area. 

Example 18. In view of equality ( 133]) it is clear that we can apply Theorem [TTJ and 
Theorem [TBI in the same way as in Example [T7J i.e. with (1 = 8 = ^ = 2 and s = 1, if the 
driving process (M t ) te [ 0T ] in ( ITS]) is a Wiener process in W 2 1 (C) with covariance operator 
Q £ i(nuc)(M^2(^))i an d even if it is a Wiener process in H\ 2 {0) with covariance operator 
Q £ £(nuc)(-f^2 2(C))- m the first case (Mt)te[o,T\ ^ oes no ^ satisfy a zero Dirichlet boundary 
condition as in Example [TT] and in the second case [M t )t£\o,T] behaves even more irregularly 
near the boundary in the sense that the first derivatives are allowed to blow up near dO. 

In these cases we choose (e K ) Ke ^ in ( 1HT|) and ( 152]) to be an orthonormal basis of the space 
W 2 (0), respectively H\ 2 {0), consisting of eigenvectors of Q G L( nnc )(W 2 (0)), respectively 
Q G £(mic)C#2 2(^))> w ith corresponding eigenvalues (A K ) Ke K. 

As in Example [13 the solution u lies in L 2 (0 T ; ^(0)) = L 2 (0 T ; -B^ 2 (0)) and, by 
Theorem [L5J it also lies in L 2 (ft T ; B£ T (0)), 1/r = a/2 + 1/2, a < 2, see Figured 

Example 19. Let the driving process (M t )te[o,T\ i n (IIS]) be a time-dependent version of 
the stochastic wavelet expansion introduced in Abramovich et al. [T] in the context of 
Bayesian nonparametric regression and generalized in Bochkina [3], ClOlCA et al. [7J- This 
noise model is formulated in terms of a wavelet basis expansion on the domain O C R d 
with random coefficients of prescribed sparsity and thus tailor-made for applying adaptive 
techniques with regard to the numerical approximation of the corresponding SPDEs. Via 
the choice of certain parameters specifying the distributions of the wavelet coefficients it 
also allows for an explicit control of the spatial Besov regularity of (M t )te[o,T] ■ We first 
describe the general noise model and then deduce a further example for the application of 
Theorem [15] 

Let {ip\ : A G V} be a multiscale Riesz basis for L 2 (0) consisting of scaling functions 
at a fixed scale level jo £ ^ an d of wavelets at level jo an d all finer levels. As in the 
introduction, the notation we use here is different from that used in Section I2T21 because we 
do not consider a basis on the whole space R d but on the bounded domain O. Information 
like scale level, spatial location and type of the wavelets or scaling functions are encoded 
in the indices A G V. We refer to Cohen Sections 2.12, 2.13 and 3.9] and Dahmen, 
Schneider [T5], [TB], [T7] for detailed descriptions of multiscale bases on bounded domains. 
Adopting the notation of Cohen we write V = Uj>j -i^j; where for j > j the set 
Vj C V contains the indices of all wavelets ip\ at scale level j and where Vj _i C V 
is the index set referring to the scaling functions at scale level jo which we denote by 
ipx, A G Vj _i, for the sake of notational simplicity. We make the following assumptions 
concerning our basis. Firstly, the cardinalities of the index sets Vj, j > jo — 1> satisfy 

C~ x 2 id < IVj-I < C2 jd , j > jo — 1. (34) 
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Secondly, we assume that the basis admits norm equivalences similar to those described in 
Theorem [3J There exists an r G IN (depending on the smoothness of the scaling functions 
ipx, A G Vj _i, and on the degree of polynomial exactness of their linear span), such that, 
given p, q > 0, max{0, d(l/p — 1)} < s < r, and a real valued distribution / G T>'(0), 
we have / G Bt, q (0) if and only if / can be represented as / = J2\ev caV-'a, (ca)agv C R 
(convergence in T>'(0)), such that 

C £ £ \c x \>) 9/P \ ^ < oo. (35) 

\i=io-i AeVj / 

Furthermore, ||/||b| (O) is equivalent to the quasi-norm f[3"5"j) . Concrete constructions of 
bases satisfying these assumptions can be found in the literature mentioned above. Con- 
cerning the family of independent standard Brownian motions (u>*) t6 [ 0) T], k e IN, in (JT]) 
respectively (fTBl) . we modify our notation and write (u^) t6 [ 0jT ], A G V, instead. The de- 
scription of the noise model involves parameters a > 0, b G [0, 1], c G R, with a + 6 > 1. 

cd a( j — ( j'q — 1) )d 

For every j > j — 1 we set a-,- = (J — (jo — 2))~2 2 and let Yx, A G Vj, be Bernoulli 
distributed random variables on (Q,J r ,'P) with parameter pj = 2~ b ^~^°~ 1 '' d , such that 
the random variables and processes Y\, (w^) t e[o,T\, A G V, are stochastically independent. 
Now we are ready to define (M t ) t6 [ ,T] by 

00 

M t := J2 a i Y ^- w t> tE[0,T\. (36) 

.5=50-1 AeVj 

Using fl35|) . fl34|) and a + 6 > 1, it is easy to check that the infinite sum converges in 
L 2 (tt T ; L 2 (0)) as well as in the space Ai^ c (L 2 (0), (Ft)) of continuous, square integrable, 
L 2 ((9)-valued martingales w.r.t. the filtration (J : 't)te[o,T\- Moreover, by the choice of the 
hyperparameters a, b and c one has an explicit control of the convergence of the infinite 
sum in (136]) in the (quasi-)Banach spaces L P2 (Qt', Bp uq ((D)), s < r, pi,q > 0, p 2 < q. 
(Compare ClOlCA et al. [7] which can easily be adapted to our setting.) 

With regard to Theorems [TTJ and [TBI let again d = p = r y = 9 = 2. Equation f)16p with 
(M t ) t€ [ ,T] defined as above corresponds to equation (jTJ if we set 

g x (u, t, ■):= o~jY\(u))ip\( • ), A G Vj, J > j - 1, (u, t) G Q T , 

and sum over all A G V instead of k G IN. In the following we write l 2 = £ 2 (V). Since 
a + b > 1 and ||#|| h ° 2 (<ot 2 ) = VVT\\ M \\L 2 (n T -,L 2 (0)) we have g G H^ 2 (£>,T;4). Let us 
impose a bit more smoothness on g and assume that a + b > 2. This is sufficient to ensure 
that g G 2 (£>, T; £ 2 ): Using flU one sees that the H^ 2 (0, T; £ 2 )-norm of # = (# A ) Ae v 
satisfies 

- 70 ngZ 
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oo 

= tv E °j Y ?WM 2 H U 0) 

j=jo-i AeVj 

oo 

L(O) 

i=io-iAeVj |a|<i 

oo 

<C E E a jPs\\MwZ(py 

Since W 2 1 (0) = B\ 2 {0) with equivalent norms we can use the equivalence (135]) with f = ip\ 
to get 

oo 

imIhi i2 (o,t ; * 2 ) < c E E^V 2 ' 

i=io-i AeVj 

oo 

= C E |V j |(j-(jo-2)) 2c 2- 2a ^-^- 1 »2- 2 ^-^- 1 »2 2j ' 
i=io-i 

oo 

<C X] (j-(jo-2)) 2c 2- 2 ^- 2 ). 
j'=io-i 

In the last step we used (1341) with d = 2. Thus (7 G H2 2 (0, T; £ 2 )- As in Example [T7] 
we may apply Theorems [TT] and [TB] to conclude that for every initial condition w G 
L 2 (f2, -^O; Pj ^22(0)) there exists a unique solution of equation (JXJ) in the class i} 2 2 ((9, T), 
which, in general, is not in /^(flr, W^iO)) for all s < 2, but it belongs to every space 
L 2 (fi T ; B rA°)) with " < 2 and r = 2/(a + 1). 

Remark 20. In practice, many adaptive wavelet-based algorithms are realized with the 
energy norm of the problem which is equivalent to a Sobolev norm. Let us denote by 
{r]\ : A G V} a wavelet Riesz basis of W^iO) for some s > 0, which can be obtained by 
rescaling the wavelet basis {ip\ : A G V} of L 2 {0), see, e.g. Cohen [8j or Dahmen [T4] , 
For the best iV-term approximation in this Sobolev norm, it is well known that 

u G B? tT {0), \ = + \ =► < CN-t*-' )/d , 

where 

0"jv,W|(O)(w) := hif - u N \\w°(G) ■ u N = ^c A r/ A : A C V, |A| < N, c x G R, A G A|. 

AeA 

Therefore, similar to the L2(0)-setting, the approximation order of the best Af-term wavelet 
scheme in W^O) depends on the Besov regularity of the object one wants to approximate. 
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wl l2 {p) 




Figure 3: Besov regularity in the scale B^ T (0), 1/r = (a — l)/2 + 1/2, vs. Sobolev 
regularity of the solution illustrated in a DeVore-Triebel diagram. 



There exist adaptive wavelet-based algorithms which are guaranteed to converge and 
which indeed asymptotically realize the convergence rate of best iV-term approximation 
with respect to the Sobolev norm. For example, Cohen, Dahmen, DeVore [5] designed 
such an adaptive numerical scheme for solving (deterministic) elliptic PDEs. First results 
for parabolic problems were obtained by Schwab, Stevenson [43J. 

Once again, the use of adaptive algorithms is justified if the rate of approximation that 
can be achieved is higher than in classical uniform schemes. Let un, N 6 IN, denote a 
uniform approximation scheme (e.g. a Galerkin approximation) of u. It is well-known that 
under certain natural conditions, see, e.g. Dahlke, Dahmen, DeVore [T2J or DeVore 
[21] or Hackbusch [26] . 

||« - u N \\ W s (0) < CN^ a - s)/d \\u\\ w?{0 ). 

This means that, even in this case, adaptivity can pay off if the Besov smoothness of the 
solution is higher than its Sobolev regularity. 

Let us discuss this relationship in more detail for the examples from above. We consider 
approximation in W^{0). As already mentioned in Example fT7l in general we cannot expect 
that the spatial Sobolev regularity of the solution is higher than 3/2. Therefore, uniform 
schemes yield an approximation rate of 0{N^ 1 ^ 4 ). 

On the other hand, our main result shows that 

u E L T (n T ; B? T (0)), - = £ + \ for all a < 2. 

Til 

Therefore, by interpolation and embedding of Besov spaces we can achieve that the solution 
is contained in all the spaces L T (Q T ; B" T {G)) corresponding to the points in the trapezoid 
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with vertices (1/2,0), (1/2,3/2), (3/2,2), (3/2,0) and to the points to the right of this 
trapezoid in the DeVore-Triebel diagram, cf. Figure [3J As a consequence, we get by a 
short computation, that 

u G L T (n T ; B« T (0)), - = ^ + \ for all a < ~. 

Thus, best iV-term wavelet approximation provides order O^N" 1 ^ 3 ), so that again the use 
of adaptivity is completely justified. 



A Convergence of the stochastic integrals 

In this section we give a proof of the M 2 /{R, (^)) -convergence of the sum of the stochas- 
tic integral processes (f (g K (s, • ), if) d> w s) te \ Ty K e ^> appearing in formula (JSJ). Let 
us assume that g G g (G , T ; £ 2 ) f° r some p G [2, oo) and 7, 9 G ft. We use an anal- 
ogous strategy to [31], Remark 3.2. Due to the independence of the Brownian motions 
(w£)te[o,T\i K G IN, the covariation process ([w K , w e ]t)t£[o,T\ vanishes if k 7^ £, and we have 
by Ito's isometry: 



E 



K = l 



/ (g K (s,-),if)dw« 



E 



E 

_K=1 

OO 

E E 

K=l 



(•) 



(g K (s, ■),<p)dw^,J2 (9 K (s, .),<p) dw K s 



K=l 



(0 



\(g K (s,-),p)\ 2 ds. 



We are going to show that the last term is less or equal a constant times \\g\ 

which is finite due to our assumption. Then the convergence of the integral processes in 
A^^ 6 (R, {J~tj) follows by Doob's maximal inequality for martingales. 

For u G V{0) and n G Z we use the notation u n := (- n {c n ■ )u{c n ■ ) G S'(R d ). 
Let us abbreviate L T (R d ) by L T for all r > 1 in the sequel. Setting p' = p/{p — 1), we 
denote by ( ■ , ■ )l p xl p , '■ L p x Zy — > R the dual form obtained by continuous extension of 
{<p, if)) = J <f(x)if>(x) dx, <f, if) G C^{K d ). Now we are ready to estimate as follows: 



00 „ 00 „ 

£ / \(g K (s,.),v)\ 2 ds = ]T / 

«=1 "n «=1 n 



ds 



T 



< 



E 

K = l 

OO 

E 

K = l 



E< 
E 



rid 1 



'nd I 



;i-A)*tf( a , 0,(1- a)- V 



L p xL pl 



ds 



A)~^ r 



1 1/2 1 



Il 2 



;i-A)-^ r 



|V2| 



I La 



c/.S' 
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oo 

K=l 



1/2 



^c 2 ^|||(l-A)?<( s ,.)|-|(l-A)-^, 



IL2 



EH(l-A)-^n|| Li ) 

\ne2 / 



1/2" 



ds. 



Here we have used Holder's inequality twice. Since tp has compact support in O and £_ n 
equals zero outside C_ n , the functions tp n vanish on R d for all but finitely many n G Z. As 
a consequence, the sum J^nez 11(1 ~~ A) -7 / 2 ^!!^ has only finitely many non-zero terms. 
Therefore, 



E f T \(g K (s,-),<p)\ 2 ds 

00 /-T 

^E / E^IKi-^feOi-Ki-A)-^ 1 / 2 !!^^ 
= c/ t Ec 2W (Eki-a)^( S ,- 

^° neZ K=l 



M(i-A)-^v 



£p/2 x -kp/(p-2) 



where the constant C depends on ip. In the last step we used the fact that 
(1 - A)^g n (s, •) = ((!- A)^ /2 ^(s, • )) KGM e L p (R d ; £ 2 ) 

P (g) A-almost everywhere in which results from g being an element of 6 (0; £ 2 ). 
Applying again Holder's inequality we obtain 



£p/2 x -kp/(p-2) 



Ec 2W ( E la - A ) 7/2 ^, • )i 2 , i(i - a)-a> b 

< e c2n i 1 (i - A ) 7/2 ^ •) 1* w'/y^ 11 (i - A)-^Vn|U p 



/(P-2) 



P-2 



<C E cn 1l(l-A) 7/2 ^, 



■)ML 



^ c 2„(d-e)p/(p-2)|| (1 _ A)-^Vn||^;^ 



<C E cn 1l(l-A) 7/2 ^, 



VngZ 



where we have used that p < 2 and that only finitely many of the ip n , n£Z, are non-zero. 
All in all we have shown 

2 

00 „ T „ T / \ p 

E/ \{g K {s,-)M\ 2 ds<C E cn H(i~ A ) 7/ V(s, OkllL) 
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Finally, taking the expectation and applying Jensen's inequality yields 



E E [ T \(9 K (s,-)^)\ 2 ds<C\\g\\^ 



(0,T;£ 2 ) 



and this finishes the proof. 



B General linear equations 

In the introduction we have indicated that our main result can be extended to equation 
02]). The major reason is, that by a result from Kim [30J an estimate similar to the one 
proved in Corollary [T5] holds not only for the model equation ([T]) but for equations of the 
type OH]), provided the coefficients a M ^,6 M , c, o^ K and r] K , the free terms / and g K and the 
initial value uq fulfil certain conditions. We can use this fact to extend our regularity result 
to such equations. In this section we want to get more precise and point out how to do 
this. 

For the convenience of the reader we begin by presenting the result from Kim [3TJ1 
Theorem 2.8]. Therefore, we need some additional notations. For x, y G O we shall write 
p(x, y) := p(x) A p(y). For a G R, 5 G (0, 1] and k G No we set: 

[f][ a) :=su V p k+a (x)\D k f(x)\ 

\ft ] ■= E[/]? } 

3=0 

whenever it makes sense. We shall use the same notations for £ 2 - va hied functions (just 
replace the absolute values in the above definitions by the £ 2 -norms). Furthermore, let's fix 
an arbitrary function 

(JL : [0, oo) -)■ [0, oo) , 

vanishing only on the set of nonnegative integers, i.e. po(j) = if and only if j G INo- We 

set 

t+ := t + Ho(t). 

Now we are able to present the assumptions on the coefficients of equation ([3]) (see Kim 
[30| Assumptions 2.5 and 2.6]). 

[Kl] For any fixed i6 0, the coefficients 

a"" (., ., x) , W (, .,x),c (., ., x) , (., ., x) , ?f (., , x) : Q x [0, T] R 

are predictable processes with respect to the given normal filtration (^t)t6[0T]- 



and [f]^ s := sup p + (x,y) — — 5 



x, y eo \x-y\ 



and := |/|^ + \ff&, 
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[K2] (Stochastic parabolicity) There are constants S ,K > 0, such that for all (cj,t,x) G 
ttx[0,T]xO and A G R d : 

5 \X\ 2 < aJ n7 (uj } t } x)X fl X u < K\X\ 2 , 

where aF 7 := - |((X M , o v )i 2 for [i,ue{l,..., d}. 

[K3] For all (w,t) G Q x [0,T]: 

K( Wl t,Olf°[ f + |^(a; I t,.)|g| + + |c(a;,t,.)lg| + 

+ K(w,t ) .)C + + Kw ) «,.)lfii| + <^ 

[K4] The coefficients a Mly and cr^ are uniformly continuous in x G O, i.e. for any e > 
there is a 5 = 5(e) > 0, such that 

|a^(w, t, x) - a^(w, t,y)\+ i, x) - ^(w, t, y) |«, < e, 

for all (a;, i) G f2 x [0, T], whenever x,y G C with |x — y| < <5. 

[K5] The behaviour of the coefficients 6 M , c and ^ can be controlled near the boundary of 
O in the following way: 

lim sup {p(x)\b >l ( y u, t,x)\ + p 2 (x)\c(uj, t,x)\ + \v(cu, t, x)\e 2 } = 0. 
pO)-k) 
xeo te[o,T] 

Here is the main result of Kim [30] . 

Theorem 21. Lei p G [2, oo) and let assumptions [Kl] - [K5] be satisfied with K,8q > 0. 
Then there exists a constant k = K (d,p, So, K, O) G (0,1) such that, if 9 G (d — Ko,d + 



ion 



K +p- 2), for any f G H^ p (C?,T), 5 G Hj/(C»,T;4) and n G U^O), e?ua* 
mt/i initial value uq admits a unique solution u G fjp (C?, T), i.e. ; t/iere exists an (up to 
indistinguishability) unique V '(C) -valued predictable process u G ETJ e (C,T), swca £/ia£ 
/or any </? G C™(0) the equality 

(u(t,.),<f) = (u(0, .),¥>) 

+ / (a^(s, .) + 6"(s, •) + c(s, >(s, .) + f(s, .), y?)ds 

oo „ t 

+ V / (o-^(s,.)u a;M (s,.) + n K (s,.)n(s,.) + ^(s,.),^)d< 
/ioWs for all t G [0, T] with probability 1. Moreover, for this solutions we have 



W<leiO,T) < C ^/ll^ (0fr) + hW^JiCTM + KHS^O) J , (37) 

where C is a constant depending only on d, 7, p, 9, 5q, K, T and the domain O . 
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An immediate consequence of this theorem is the following estimate. 

Corollary 22. In the situation of Theorem{21\ with 7 = m G IN, the following inequality 
holds for every r G [0,p]: 

p m -P\D m u(u,t, ■)W\\l P {0)dtnd^) 

< C (ll/lle™- + 2 p + hWn^iCT-M + \\ u o\\u™ (o)) , 

where = 1 + 

Proof. Just repeat the arguments of the proof of Corollary [121 and use estimate (1371) instead 
of (JToT) at the beginning. □ 

Now we can present our main result in the generalized setting. 

Theorem 23. Let p G [2, 00) and let assumptions [Kl] - [K5] be satisfied with appropriate 
constants K,5 > 0. Moreover, let f G lT p ~l p (0,T), g G W~g{0,T\l 2 ) and u G U2 fi {0) 
for some 7 G IN. Denote by u the unique solution of equation (J3J) in the class Sj pd (0,T) 
for a given 6 G (d — Ko,d — 2 + p + Kq), where k q = K (d,p, 80, K, O) G (0, 1) as in Theorem 
IM Assume furthermore that u G L p (Q T ; B s pp (0)) for some s G (0,7 A (1 + ^)]. 
Then, we have 

ueL T ({l T ;B« T (0)), - = ^ + -, for all «g(o, 7 A-^-), 

rap \ a — 1/ 

and the following estimation holds 

|MU T (n T ;B? iT (0)) < C (H/HhJ-^^t) + II5 , Hh^- 1 (C>,T;£ 2 ) + 1 1 U 1 1 U^ g (O) + 1 1 ^ I U p (n r ;J B|, p (C>)) ) • 



Proof. We can argue like we did in the proof of Theorem [T51 We just have to use Corollary 
where we used Corollary [T2J □ 
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